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ABSTRACT 

We report the detection of the polar CE Gru in X-rays for the first time. We 
find evidence for a dip seen in the hard X-ray light curve which we attribute to the 
accretion stream obscuring the accretion region in the lower hemisphere of the white 
dwarf. The X-ray spectrum can be fitted using only a shock model: there is no distinct 
soft X-ray component. We suggest that this is because the reprocessed component is 
cool enough so that it is shifted into the UV. We determine a mass for the white dwarf 
of ~1.OM . 
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1 INTRODUCTION 

CE Gru (also known as Grus V-l and Hawkins V-l) was 
discovered by Hawkins (1981, 1983) during a search for vari- 
able objects using UK Schmidt plates. Further observations 
, by Tuohy et al (1988) confirmed its binary nature (an or- 
bital period of 108.5 mins) while Cropper et al (1990) found 
strong circular polarisation - a characteristic of the polar (or 
AM Her) class of cataclysmic variable. These objects are in- 
teracting binaries in which material flows from a dwarf main 
sequence star onto a magnetic (B ~10-200MG) white dwarf. 
This strong magnetic field is high enough to force the spin of 
the white dwarf to be synchronised with the binary orbital 
period. 

Tuohy et al (1988) demonstrated that there are two ac- 
cretion poles which are characterised by different emission 
properties: one pole which is always in view and is stronger 
in blue light, while the other pole which is visible for only 
~0.35 of the white dwarf spin period is stronger in red light. 
Cropper et al (1990) showed that the pole always in view 
was positively circularly polarised, with the other being neg- 
atively polarised. 

The UK Schmidt plates show that CE Gru exhibits two 
distinct levels of brightness, (B ~18 and B ~21). In the 
fainter state the accretion flow is much reduced, or stopped 
altogether. This may be the reason why CE Gru was not 
detected during the ROSAT all-sky survey (Verbunt et al 
1997). Indeed, CE Gru has never before been detected in 
X-rays. In this paper, we report the first X-ray detection of 
CE Gru which were made as part of a survey of polars using 
XMM-Newton. 



2 OBSERVATIONS 

The satellite XMM-Newton was launched in Dec 1999 by 
the European Space Agency. It has the largest effective area 
of any X-ray satellite and also has a 30 cm optical/UV tele- 
scope (the Optical Monitor, OM: Mason et al 2001) allowing 
simultaneous X-ray and optical/UV coverage. CE Gru was 
observed using XMM-Newton on 2001 Oct 31. The EPIC 
instruments (imaging detectors covering the energy range 
0.1-10keV with moderate spectra resolution) were operated 
in full frame mode (the count rate was not high enough to 
cause pile-up problems). The RGS detectors (high resolu- 
tion spectrographs operating in the 0.3-2.0keV range: den 
Herder et al 2001) were configured in the standard spec- 
troscopy mode. We clearly detect CE Gru in the X-ray band. 
OM data were taken in two UV filters (UVW1: 2400-3400 
A, UVW2: 1800-2400 A) and one optical band (V band). 
The observation log is shown in Table [l]. 

The data were processed using the XMM-Newton Sci- 
ence Analysis Software v5.2. The RGS spectra were of low 
signal to noise and showed no evidence for significant line 
emission: we do not consider them further. For the EPIC 
pn detector (Striider et al 2001), data were extracted us- 
ing an aperture of 40 arc sec centered on the source. 
Background data were extracted from a source free region. 
For the EPIC MOS detectors (Turner et al 2001) we ex- 
tracted data in a similar way, but extracted the background 
from an annulus around the source. The background data 
were scaled and subtracted from the source data. In ex- 
tracting the EPIC pn spectrum, we used only single pixel 
events and used the response file epn JT20_sY9_thin.rmf. In 
the case of the MOS data we used the response files m[l- 
2]_thinlv9ql9t5r5_all_15.rsp. The OM data were analysed in 
a similar way using omichain and omf chain (this latter task 
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Detector 


Mode 


Filter 


Exp (s) 


EPIC MOS 


Full frame 


thin 


7792 


EPIC pn 


Full frame 


thin 


5468 


RGS 


Spectroscopy 




8255 


OM 


Image / fast 


UVW1 


1500 


OM 


Image / fast 


V 


1500 


OM 


Image / fast 


UVW2 


1800 



Table 1. The log of XMM-Newton observations of CE Gru. 

was not incorporated in SAS v5.2 but will be in a later ver- 
sion). Data were background subtracted and corrected for 
coincidence losses (Mason et al 2001). 

The Optical Monitor data shows that CE Gru had a 
mean brightness of V=17.9 and a maximum V=17.5. Since 
the V band observations started at the descent from max- 
imum its likely that its true maximum was brighter than 
this. Assuming a similar colour to that found by Tuohy et 
al (1988) (B - V=0.53) this places CE Gru in a high ac- 
cretion state at the time of the XMM-Newton observations 
and a similar brightness to observed by Tuohy et al (V ~18). 
The length of the observation in the EPIC MOS detectors 
covered just over 1 orbital cycle. 

The flux in the UV filters corresponds to: UVW1 
~ 4.5 x 10~ 16 ergs s" 1 cm- 2 A, UVW2 ~ 4.0 x 10" 16 
ergs s _1 cm _2 A, (based on OM observations of isolated 
white dwarfs). The mean flux in the V filter corresponds 
to ~ 2.5 x 10" 16 ergs s" 1 cm" 2 A at 5000 A. 



3 LIGHT CURVES 

We show in Figure [j] the light curve of CE Gru in various 
energy bands using EPIC pn and EPIC MOS data (where 
the MOS1 and MOS2 data have been co-added) folded on 
the orbital period of Tuohy et al (1988). We also show the 
OM data: CE Gru was detected in all three filters (it is also 
the first time that it has been detected in the near-UV). 

There are two distinct parts to the orbital light curve: a 
faint phase lasting ~ 0.6 orbital cycles, and a brighter phase 
lasting ~ 0.4 orbital cycles. These relative durations are sim- 
ilar to those in the optical light curves seen by Cropper et al 
(1990). The bright phase is also of a similar duration to that 
of the red pole seen by Tuohy et al (1988) when the system 
was at a similar brightness compared to our observations. 
We therefore assign the bright phase X-ray emission to the 
pole in the lower hemisphere, and the fainter phase emission 
to the pole in the upper hemisphere. Although there maybe 
some ambiguity in this assignment, the fact the we observe 
decreasing emission in the V band at the end of the bright 
phase provides supporting evidence for this assignment since 
the 'blue' pole does not show such a rapid drop in V. 

In the bright phase, a dip is seen in the soft X-ray light 
curve, but not at higher energies. This is characteristic of 
photo-electric absorption, and is seen in other polars (eg 
Watson et al 1989). It is thought to result from the accre- 
tion stream crossing our line of sight to the accretion region, 
thereby absorbing the soft X-rays. This is consistent with 
the model of Wickramasinghe et al (1991) for CE Gru in 
which the two accretion regions lie near the magnetic poles, 
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Figure 1. The light curves obtained of CE Gru using XMM- 
Newton. The X-ray data were obtained using the EPIC pn and 
MOS detectors. The MOS1 and MOS2 light curves have been 
co-added. The Optical Monitor data (bottom plot) shows the op- 
tical/UV data with units of flux of 10 -16 ergs s" 1 cm" 2 . The data 
have been phased on the orbital period of Tuohy et al (1988) and 
binned into S<j> =0.02 bins. We have chosen phase 0.0 to corre- 
spond to the dip centered in the bright phase. 



approximately at the foot-points of the field lines passing 
through the region where the stream threads onto the mag- 
netic field. It is most likely that the dip is caused by the 
accretion stream to the upper pole absorbing the emission 
from the lower pole: there is clear evidence for accretion at 
the upper pole and it is inevitable that the stream crosses 
our line of sight. It is possible for the ballistic stream or the 
stream to the lower pole to be the cause of the absorption, 
but this requires a high inclination and the ballistic stream 
to penetrate very close to the white dwarf before threading. 
These are special conditions, and in the absence of even a 
grazing eclipse we consider this to be unlikely. 

At energies greater than 2keV, there is no evidence for 
the dip seen at lower energies. Using the spectral model de- 
scribed below we can estimate the absorbing column density 
of the accretion stream which causes the dip by increasing 
the column until we match the observed count rate at lower 
energies. We estimate that the total column density in our 
line of sight must be ~ 5 — 10 x 10 21 cm -2 . This is of the 
same order as that seen in other polars with stream dips 
(Watson et al 1989). 



4 X-RAY SPECTRA 
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4.1 The model 

We extracted a faint and bright phase spectrum (0=0.2-0.8 
and 0=0.8-1.2 respectively) from the EPIC MOS data and 
a faint phase spectrum from the EPIC pn data (these data 
covered too short a length of time to obtain a useful bright 
phase spectrum). We modelled the data using a simple neu- > 
tral absorber and an emission model of the kind described < 
by Cropper et al (1999). This emission model, unlike single .» 
temperature thermal bremsstrahlung models, is a realistic | 
physical description of the post-shock accretion region in ° 
polars. It is based on the prescription of Aizu (1973) which J 
predicts the temperature and density profiles over the height e 
of the accretion shock. However, it has been modified to take 
into account cyclotron cooling (which can be significant in 
polars) and also the variation in gravitational force over the 
shock height. To reduce the number of free parameters we 
fix the parameter, e s , (the ratio of cyclotron cooling to ther- 
mal bremsstrahlung cooling), at 5: this implies a magnetic 
field strength typical of polars (20-50MG). We also fix the 
specific accretion rate at 5 g s _1 cm -2 (typical of polars in 
a high accretion state). Changing these parameters does not 
have a great effect on the results for data with low-moderate 
signal to noise ratio (as in our CE Gru data). Observations 
of polars using ROSAT showed that many polars had promi- 
nent blackbody components (Ramsay et al 1994). Therefore, 
we also added a cool blackbody model and determined the 
change in the fit. 

4.2 Results 

We initially consider the EPIC spectra from the faint state. 
We show in Table ^| the fits to the MOS and pn spectra (we 
fit the MOS spectra simultaneously while leaving their nor- 
malisation parameters untied). We find that the spectrum 
can be well modelled without a blackbody component. The 
Hydrogen column density is low (< 10 20 cm -2 ). Adding a 
blackbody to the model makes no significant difference to 
the fit. We do, however, show in Table H the fit and flux when 
we add a blackbody of temperature 30eV (typical of the 
temperature of the blackbody component seen in ROSAT 
observations of polars) to the model. In the bright phase 
data, again we find no significant improvement to the fit 
when we add a blackbody to the model. We show in Figure 
^| the EPIC MOS1 spectra taken from the bright and faint 
phase, together with the best fits (assuming no blackbody). 

Table ^ shows that the inferred mass of the white 
dwarf (assuming the white dwarf mass-radius relationship 
of Nauenberg 1972) is ~1.0M Q . Although this is higher 
than found for isolated white dwarfs, it is comparable with 
the mass of other accreting magnetic white dwarfs. Ramsay 
(2000) show that the latter class are biased towards higher 
masses. 

4.3 Soft and hard luminosities 

We define the hard X-ray luminosity as (Lhard = 
47rFlux, lt „. di f, id 2 ) where FIuxw^m is the unabsorbed, 
bolometric flux from the hard component and d is the dis- 
tance. Since a fraction of this flux is directed towards the 
observer, we switch the reflected component to zero after 
the final fit to determine the intrinsic flux from the optically 
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Figure 2. The EPIC MOS1 spectra extracted from the bright 
phase (dots) and the faint phase (diamonds) data together with 
a multi-temperature shock model (solid line for the bright phase 
and dashed line for the faint phase). 



thin post-shock region. We define the soft X-ray luminosity 
as (L so ft = 7rFlux S0 /t,b ;sec((9)d 2 ), where we assume that 
the soft X-ray emission is optically thick and can be ap- 
proximated by a small thin slab of material, the unabsorbed 
bolometric flux is Flux so j- ti (, ; and 8 is the mean viewing 
angle to the the accretion region. 

The is some degree of uncertainty in the viewing angle 
to the accretion regions. Tuohy et al (1988) found i ~ 40° 
and the red pole (the bright X-ray pole) has /3 ~ 125° and 
the blue pole (the faint X-ray pole), ft ~ 40°. More detailed 
modelling of Cropper et al (1990) and Wickramasinghe et al 
(1991) find i ~ 40°, a dipole offset of 165° and the red pole 
displaced from the magnetic pole by ~ 30° and the blue pole 
~ 150°. However, the blue accretion region (the faint X-ray 
region) is seen for all spin phases (and hence a relative small 
viewing angle) and the red region for only a small fraction 
of the spin phase (and hence at a high viewing angle). For 
argument we apply a mean viewing angle of 0=60° (implying 
sec#=2.0) for the faint phase and #=80° (sec#=5.8) for the 
bright phase. 

We show in Table ^| the resulting luminosities for the 
shocked component (the 'hard' component) and the black- 
body component assuming fcT(,j,=30eV. We find that the ra- 
tio, L so ft/ Lhard, in the faint phase is very low, with an up- 
per limit of ~0.1. In the bright phase, the upper limit is 
0.7. In the standard accretion model, hard X-rays in the 
post-shock region irradiate the photosphere of the white 
dwarf and are re-emitted at lower energies. Assuming the 
re-processed component is emitted in soft X-rays, we expect 
Lsoft/ L s hock ~ 0.5 (Lamb & Masters 1979, King & Lasota 
1979). The total shock luminosity is the sum of hard X-ray 
emission and the emission from the cyclotron component 
which originates from the post-shock flow as well. There- 
fore, L h ard underestimates L shock- 

In the faint phase, it is clear that for a blackbody of 
temperature 30eV, the resulting energy balance is not con- 
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Faint Phase 


Bright Phase 






EPIC MOS 


EPIC pn 


EPIC MOS 


N H (10 20 cm- 2 ) 




0.0 + J 


0.0 


00 +0.02 


Afl (Mr-,') 




L04 -0.23 


1 03+° 18 


1 29 n 1 s 
J--*<E?— 0. 15 




r cm 2 ) 


o i +0 9 i n-12 
2-l_ .6 x iU 


9 9 +0.5 -12 
z - z -0.3 x 


7 l+°' 8 x 10" 12 






2.1+1;° x 10 -i2 


8.2+ ; 2 x 10- 12 


L h ard,bol (ergs S _1 


^100 ) 


9 c^+l- 1 v 1 fl^O 

Z.O q y X _1_U 

2.6lo' 2 x 10 30 


9 7+O.5 v -1 ri30 
z - ' —0.4 x iu 


9.lf°J x 10 30 
I.oIq 2 x 10 31 


Xg(dof) 




1.07 (82) 


1.08 (104) 


0^92 (143) 


+ bb (eV) 




30 


30 


30 


Flux so/t:bo ( (ergs s 


- 1 cm -2 ) 


1.5+H x 10" 13 


1.9li 7 g x 10" 13 


l.l+ 2 '5 x 10- 12 






1.5+j' 2 x 10~ 13 


8.9±i 4 9 6 x 10- 13 


L ao ft,bol (ergs s" 1 


^100 ) 


9.0+H x 10 28 
9.0±i; 8 x 10 28 


l.ltH x 10 29 


1.9+fg x 10 30 
x !0 30 


xS(dof) 




1.06 (79) 


1.07 (103) 


0^90 (140) 



Table 2. The parameters for the EPIC spectra extracted from the faint phase and bright phase. The MOS1 and MOS2 spectra were 
fitted simultaneously with all the parameters linked except the normalisations. We also show the effect on the fit if we add a blackbody 
with fcT(,i,=30eV to the model. 



k l'iii xt Blackbody flux Blackbody L so f t J 

(eV) (dof) (ergs s -1 cm -2 ) luminosity Lhard 

(ergs s" 1 d\ m ) 

20 4.31 (27) 2.68xl0~ 10 1.6xl0 32 64 

10 1.39 (27) 3.86X10- 11 2.2xl0 31 9.0 

5 1.43 (27) 6.70X10™ 12 4.0xl0 30 1.6 

2 1.44 (27) 1.34X10" 12 8.0xl0 29 0.3 



Table 3. Using the MOS spectra from the faint phase we have 
added a blackbody of various temperatures: the normalisation 
has been set to give the observed flux in the UVW2 filter. We 
show the resulting fit to the X-ray data in the 0.1-1.0keV band, 
the unabsorbed bolometric flux and luminosity of the blackbody 
component and the ratio L so f±/ Lhard (cf Table ^). 

sistent with the standard shock model. However, lower tem- 
perature blackbodies can be 'hidden' at EUV wavelengths. 
To investigate this further, we added a blackbody compo- 
nent of different temperatures and set their normalisation so 
that the derived flux matched that observed in the UVW2 
filter. We show the model of the combined blackbody plus 
post-shock model in Figure ^ for fcTt,{,=2, 5, 10 and 20eV. 
We also show in Table tfl the fit to the data for these models 
in the 0.1-1.0keV band, the resulting flux and luminosity 
for the unabsorbed blackbody. We find that the observed X- 
ray spectra are not consistent with blackbody temperatures 
greater than ~10eV. 

The un-heated surface of the white dwarf is also ex- 
pected to contribute to the flux seen in the UV. However, 
we do not know its temperature (but it is expected to be 
near l-2eV ~10000-20000K) nor do we have an accurate es- 
timate of its distance. Because of these uncertainties we do 
not add a blackbody to account for this component: ~10eV 
is therefore an approximate upper-limit to the temperature 
of the reprocessed component. For fcTi,6=10eV the result- 
ing L so ft /Lhard ratio (~9) is very much greater than that 
expected from the standard accretion shock model but for 
fcT;,(,=5eV it is slightly higher than expected (~1.6). For 
lower temperatures (2eV) the ratio is slightly lower than 
that expected. 



4.4 The distance 

In the bright phase the hard X-ray bolometric luminosity 
is ~ 1 X 10 31 ergs s _1 dioo- To place a very crude estimate 
on the distance to CE Gru we compare this luminosity to 
the hard X-ray luminosity of AM Her (whose distance is 
reasonably well determined). Ishida et al (1997) observed 
AM Her using ASCA when it was in a high accretion state 
and find Lhard = 1.6 x 10 32 ergs s _1 using a distance of 75 
pc. We find that CE Gru would have to lie ~400 pc distance 
to equal the Lhard found for AM Her. 



5 DISCUSSION 

Although CE Gru was detected in X-rays for the first time 
and was clearly in a high accretion state, no distinct soft 
X-ray component was observed. On the face of it, this is 
surprising since a strong soft X-ray flux has long been con- 
sidered one of the defining properties of polars. Observations 
of polars using EXOSAT (eg Osborne 1988) and ROSAT (eg 
Ramsay et al 1994) found a strong distinct soft X-ray com- 
ponent, and in general, the ratio L so f t / Lhard 3> 1. Further, 
the ratio was correlated with the magnetic field strength of 
the white dwarf: a high field gave a high ratio. To account 
for the large 'soft X-ray excess' seen in many polars, dense 
'blobs' of material which do not form a shock and radiate in 
soft X-rays were proposed (Kuijpers & Pringle 1982, Frank, 
King & Lasota 1988). 

Currently, four polars have been observed using XMM- 
Newton. Of those systems observed in an intermediate or 
high accretion state, WW Hor showed no distinct soft X-ray 
component (Ramsay et al 2001) while BY Cam has one pole 
which did show such a component and one which did not 
(Ramsay & Cropper 2002). Indeed only one other system, 
DP Leo, has shown a distinct soft X-ray component (Ramsay 
et al 2001). Even for those polars which showed the lowest 
ratios using ROSAT (EF Eri and AM Her) a distinct soft 
X-ray component was still observed. 

A variation of the standard shock model has been pro- 
posed by Heise & Verbunt (1988) and Gansicke, Beuermann 
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Energy Energy 

Figure 3. We show the best fit shock model that was used to fit the EPIC MOS1 data in the faint phase together with a blackbody 
of various temperatures. The normalisation of the blackbody was chosen to fit the observed flux in the UVW2 filter. A reprocessed 
component with less than kT^ ~10eV is therefore be cool enough not to be detected in the X-ray band. 



& de Martino (1995) who suggest that the reprocessed X- 
ray component lies in the EUV band while the strong soft 
X-ray component seen in many polars originates from dense 
'blobs' of material which do not form an exposed shock. Our 
data are consistent with this view if the temperature of the 
reprocessed component is ~ 2 — 5eV: below 2eV the flux of 
the reprocessed component is too small, while above 5eV it 
would be too high. A component with a temperature above 
lOeV would be evident in the spectrum (Figure^). The im- 
plication in this case is that the fraction of blobby accretion 
is small in CE Gru. Again this is consistent with the un- 
folded light curves (not shown) which do not show strong 
flaring seen in many systems (such as BY Cam - Ramsay & 
Cropper 2002). 

The amount of irradiation that the white dwarf receives 
is a function of height in the post-shock region, with greater 
temperatures from higher up, but more overall flux (even at 
higher energies) from near the base (Cropper, Wu & Ram- 
say 2000). On the other hand, because of the curved sur- 
face of the white dwarf, the illumination is decreasing, with 
slightly more than the square of the distance from the axis 
of the post-shock region. Further, albedo varies as a func- 
tion of temperature. Although the effects of irradiating a 
white dwarf atmosphere have been explored to some extent 



(Williams et al 1987, Heise 1995) , this is an area which needs 
further work to determine how the temperature of the re- 
processed spectrum is effected by parameters such as the 
accretion rate and magnetic field. 

We note, however, that higher specific accretion rates 
result in lower shock heights, and higher magnetic fields also 
reduce the height of the post-shock region (see for example 
Cropper et al 1999), so the solid angle of emission from the 
post-shock flow intercepted by the white dwarf photosphere 
is increased. In these high-state data of CE Gru, we have no 
reason to expect anomalously low specific accretion rates, 
and CE Gru has a magnetic field strength typical of polars. 
Therefore there is no obvious reason for the reprocessed com- 
ponent to have moved into the EUV if it is normally in the 
soft X-ray band in polars. This gives support to the Heise & 
Verbunt (1988) suggestion that the soft X-ray component is 
indeed caused by blobby accretion. 

Some factor(s) must determine the number, length and 
density of blobs. Obvious parameters include the magnetic 
field strength and orientation of the white dwarf, the or- 
bital period and the mass transfer rate. Ramsay et al (1994) 
identified the magnetic field strength as one important pa- 
rameter. Further progress in this regard awaits a systematic 
analysis of the strength of the soft X-ray component as a 
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function of these parameters in a sufficiently large sample of 
polars. 
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